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Abstract. Strongly nonlinear phononic crystals were assembled from stainless steel spheres. Single solitary waves and 
splitting of an initial pulse into a train of solitary waves were investigated in different viscous media using motor oil and 
non-aqueous glycerol to introduce a controlled viscous dissipation. Experimental results indicate that the presence of a 
viscous fluid dramatically altered the splitting of the initial pulse into a train of solitary waves. Numerical simulations 
qualitatively describe the observed phenomena only when a dissipative term based on the relative velocity between particles 
is introduced. 



INTRODUCTION 

In recent years, many experimental and numerical ef- 
forts have been focusing on the propagation of strongly 
nonlinear solitary waves in granular media [1]. These 
waves are a natural extension of the well known weakly 
nonlinear solitary waves such as the Korteweg-de Vries 
(KdV) solitons. Chains composed of steel, brass, glass, 
Nylon [1, 2, 3] and Polytetraflouroethylene (PTFE) [4] 
particles support solitary waves using Hertz' law for par- 
ticle interactions. 

However, dissipation plays a significant role in the 
transmission of pulses in almost all experimental set- 
tings. Restitution coefficients and velocity dependent 
friction have been used to investigate dissipation in a 
chain of particles [5, 6]. The energy losses are signifi- 
cant even when the experiments are performed in air (see 
Fig. 1.25 in [1]). These losses may be attributed to un- 
controlled features in the experimental setup as well as 
inherent material properties [7]. 

To predictably control the dissipation, a chain of 
spheres may be immersed in viscous fluids. The viscous 
dissipation of pulses in chains of particles differs from 
a two particle interaction in liquid because the compres- 
sion wave dominates the system's dynamic behavior. To 
our knowledge, there has not been an attempt to extend 
the research involving the collision of particles in fluids 



to a chain of particles. This paper extends current knowl- 
edge of a two particle collision in a viscous medium to a 
case where multiple particle interactions support a soli- 
tary wave. 



ANALYTICAL MODEL AND 
NUMERICAL CALCULATIONS 

The analytical model for a one-dimensional Hertzian 
chain of spheres can be found in [1]. The description of 
particle interactions in fluids is presented in [8-12]. One 
of the complications of using current sphere-fluid models 
for a chain of spheres is that a developed flow around the 
sphere is assumed before the particle collision though 
the duration of the pulse in a chain is relatively short 
(~50 /is). Nevertheless we use the outlined approach 
for the first stage of our research. In [10] two particles 
are considered to be traveling towards each other in a 
surrounding medium. If the spheres are the same size 
then the total kinetic energy is 
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where m = ^npP? is the mass of the particle, m' = 
jjTipfR^ is the added mass, R is the particle's radius, h is 



the separation distance between particle centers and U — 
is the relative particle velocity. Differentiating Eq. (1) 
with respect to h yields the added mass and pressure 
force terms. The equations of motion for a chain of 
spheres, with mass ra,-, placed vertically in a gravitational 
field in a fluid becomes 
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where = (Ri +Ri+i +Xj — jc,-+i) and the x,'s are the 
particle positions. The compression force F c j is based on 
Hertz law between particle '/' and the adjacent particles 
'z-1' and'j+1'; 

F c ,i = ^(8i. hi )-Y(8 i!i+l ), (3) 

where and \j/ can be expressed as 
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and 
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The coefficients in Eqs. (4) and (5) are identical except 
for a shift of indices. The equation for A,_i , can be 
inferred from 
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where E, and v, are the elastic modulus and Poisson's 
ratio of the particle. For all of the calculations performed 
with air as the surrounding medium, only the first two 
terms on the right hand side and the first term on the left 
side of Eq. (2) are used. When liquids surround the chain 
of spheres the buoyancy, drag, pressure, added mass and 
dissipative terms are used. The buoyancy force for each 
particle is 

F bt i = --7tR?p f g, (7) 

where g is the gravitational constant, is the radius of 
the particle and P/ is the density of the fluid. The drag 
force has a correction factor to account for a Reynold's 
number Re greater than unity i.e. not in the Stoke's flow 
regime, 

F DJ = -6%vRiUi(\ +0.15fle ' 687 ), (8) 

where v is the dynamic viscosity and Ui is the particle 
velocity. The pressure force in [9] is written for one 
sphere moving toward a stationary sphere. In our case, 
we use a relative velocity between particles, 

F pd = \nR 2 p f [(Ui-i - Ui) 2 - (U, - U i+l f] . (9) 
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FIGURE 1. Single solitary wave in calculations (a) and ex- 
periments (c) in a chain surrounded by air. Results for an iden- 
tical chain surrounded by glycerol in calculations (b) and ex- 
periments (d). Vertical scale is 2 N/div. 



We introduce an additional dissipative term based on 
the relative velocity between particles with a fitting pa- 
rameter c due to the lack of a qualitative agreement 
between the experiments and calculations based on the 
drag force term (Eq. (8)). This addition is tantamount to 
adding a dash-pot between neighboring particles [13], 

F d . i = c(U i - 1 -2U i + U i+l ), (10) 

where the coefficient c is a fitting parameter. The physi- 
cal reason for this term can be due to the radial flow of 
liquid caused by the change of contact area between par- 
ticles. 

MATLAB's intrinsic ODE45 solver was used to march 
the explicit calculation forward in time with a time-step 
of 0.05/xs. The error in the energy calculations were 
found to be less than 10~ 9 % in air and within 10~ 5 % in 
the fluid. The error in the conservation of linear momen- 
tum performed with a chain in air was less than 1(T 12 %. 



EXPERIMENTAL PROCEDURES AND 
RESULTS 

In the experiments, the impulse propagation was in- 
vestigated in three different media: air, SAE 10W-30 mo- 
tor oil, and non-aqueous Glycerol GX0185-5. The den- 
sity and dynamic viscosity used in numerical calcula- 
tions were pf = 880 kg/m 3 , v = 0.067 Ns/m 2 for oil 
and pf = 1260 kg/m 3 , v = 0.62 Ns/m 2 for glycerol. An 
experiment was performed with each of the three types 



of media to see how single and multiple solitary waves 
propagate through the chain. The chain was placed into 
an adjustable holder that had four contact points on each 
sphere. The air or fluid was able to flow freely between 
the contacts as opposed to the cylindrical holder that had 
been used in previous experiments. To create a single 
solitary wave, a spherical stainless steel striker with a 
radius of R = 4.76 mm and a mass of m = 0.4501 g 
was used to impact the top of a chain of 19 stainless 
steel particles (also with R = 4.76 mm) with a velocity of 
U() = 0.44 m/s. To create multiple solitary waves in the 
same chain, a cylindrical alumina striker with a larger 
mass of m = 1.23 g impacted the chain at Uq = 0.44 
m/s. The elastic modulus and Poisson's ratio of the stain- 
less steel particles were 193 MPa and 0.3, respectively. 
The experimental results were recorded via piezoelectric 
gauges imbedded [4] in the 10th and 15th particles from 
the top of the chain. The procedures outlined in [4] were 
implemented to compare the calculations to the experi- 
mental results using the averaged dynamic force as noted 
on the vertical axis of Fig. 1, 2 and 3. 

In Fig. 1 the numerical and experimental results are 
shown for a single solitary wave in a chain surrounded 
by air (a), (c) and glycerol (b), (d). It is evident that there 
is a very small difference between the numerical results 
for the chain in air and glycerol using Eqs. (2)-(9), which 
prompted the inclusion of the additional dissipative term 
(Eq. (10)). Without this dissipative term, the numerical 
results for glycerol shown in Fig. 1(b) do not exhibit the 
shape and speed of the pulses in experiments. 

In experiments the speeds of the single pulse were 
V s = 520 m/s in air and increased to V s = 541 m/s in 
glycerol. It is interesting to note that the calculated signal 
speed was V s = 564 m/s in air and decreased to V s = 540 
m/s in glycerol. In numerical calculations the speed of 
the single pulses in chains surrounded by both air and 
glycerol should be higher than the experimental pulse 
speeds due to the higher amplitude of the waves. Also, 
one would intuitively think that the pulse speed in air 
would be higher than in glycerol due to viscous dissi- 
pation. This is the case in calculations (even when the 
additional dissipative term Fj is added) but the opposite 
is true in experiments (Fig. 1(c) and (d)). 

In Fig. 2 the numerical results for a train of solitary 
waves in air and glycerol are presented. The dissipative 
term F c / was not included in the calculations in (a) and 

(b) and it is apparent that the amplitudes of the waves 
are much higher than the experimental results in glycerol 

(c) . Additionally, no signal splitting into a train of soli- 
tary waves were present in experiments Fig. 1(c), which 
is not reflected by the calculations without the additional 
dissipative term Fj. Again, the wave speeds in experi- 
ments were higher in glycerol than in air. 

There is a disparity between the presented analytical 
formulation using the corrected Stokes drag for dissi- 
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FIGURE 2. (a) Numerical results for multiple solitary waves 
in a chain of 19 particles surrounded by air. (b) Numerical 
results for multiple solitary waves in an identical chain sur- 
rounded by glycerol, (c) Experimental results related to (b). 
Vertical scale is 2 N/div. 
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FIGURE 3. (a) Numerical results for a single pulse with 
additional dissipative term Eq. (10) in a chain surrounded by 
glycerol, (b) Numerical result for an impact by an alumina 
striker on the same chain. Vertical scale is 2 N/div. 



pation and the experiments. The largest dissipative term 
in the calculations (before the addition of the dissipative 
term F £ /) is the drag force term. The reduction of the am- 
plitude of a single pulse was about 3% in oil and about 
4% in glycerol when using Eqs. (2)-(9). These equations, 
while important for accurately describing a particle tra- 
jectory pre and post collision are negligible when then- 
effect on the compression wave are examined. 

The results of adding the dissipative term Eq. (10) are 
presented in Fig. 3 for single and multiple pulses travel- 
ing in a chain submersed in glycerol. A similar asymme- 
try and widening of the pulse, apparent in experiments, 
can also be seen when Fig. 3(a) is compared to Fig. 1(d). 
The lower line in Fig. 1(d) has a noticeable shock-like 
tail. In calculations the tail of the wave becomes more 
shock-like as parameter c increases. The qualitative be- 
havior of the pulse in experiments is reproduced in calcu- 
lations by adding an additional dissipative term F c / with a 
coefficient c = 6.0 Ns/m for glycerol and c = 0.648 Ns/m 
for oil. These values provided the best comparison be- 
tween experimental and numerical data and were scaled 
according to the difference in viscosity. When comparing 
Fig. 2(b) and Fig. 3(b), notice the additional dissipation 



Eq. (10) has also created a shock-like response of the 
train of solitary waves. In calculations the amplitude and 
tendency to split decreased by adding this term in accord 
with experiments. 



CONCLUSIONS 

The experimental results indicate a qualitative change 
of the propagating shock and solitary waves in a chain 
immersed in glycerol while only a small change in oil. 
Without the relative velocity based dissipative term Fj, 
the equations pertaining to the surrounding fluid could 
not accurately reproduce the amplitude or the shock like 
structure of the incident pulse. This term provided the 
qualitative change needed to match the numerical analy- 
sis to the experiments. The numerical analysis predicted 
a decrease in solitary wave speed as the viscosity of 
the surrounding fluid increased contrary to experiments. 
This phenomenon may be explained by an increased ef- 
fective stiffness modulus between particles in the pres- 
ence of a viscous fluid. 
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